Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 10 September 2012 (MN WTp^ style file v2.2) 



The properties of the star-forming interstellar medium at 
z = 0.84—2.23 from HiZELS: Mapping the internal dynamics 
and metallicity gradients in high-redshift disk galaxies 



10 September 2012 



^:A.M. Swinbank,!'* D. Sobral,^ Ian Smail,i J. E. Geach,^ P. N. Best,^ I. G 
ClI McCarthy,^ R. A Crain,^ & T. Theuns^'^ 

^ Institute for Computational Cosmology, Durham University, South Road, Durham, DHl 3LE, UK 
^Leiden Observatory, Leiden University, PC Box 9513, 2300 RA Leiden, the Netherlands 
' ^Department of Physics, McGill University, Ernest Rutherford Building, 3600 Rue University, Montreal, Quebec H3A 2T8, Canada 
'^SUPA, Institute for Astronomy, University of Edinburgh, Edinburgh, EH19 3HJ, UK 
I I, ^School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, UK 
^'"'j ■ ^University of Antwerp, Campus Groenenborger, Groenenborgerlaan 171, B-2020 Antwerp, Belgium 
^ ' "Email: a.m.swinbank@durham.ac.uk 

^' 

6 

^: 

. ABSTRACT 

We present adaptive optics assisted, spatially resolved spectroscopy of a sample of nine 
Ha-selected galaxies at = 0.84-2.23 drawn from the HiZELS narrow-band survey. 
These galaxies have star-formation rates of l-27M0yr~^ and are therefore represen- 
tative of the typical high-redshift star-forming population. Our '--^kpc-scale resolution 
I observations show that approximately half of the sample have dynamics suggesting 

■ that the ionised gas is in large, rotating disks. We model their velocity fields to infer 

the inclination-corrected, asymptotic rotational velocities. We use the absolute i?-band 
magnitudes and stellar masses to investigate the evolution of the _B-band and stellar 
mass Tully-Fisher relationships. By combining our sample with a number of similar 
I measurements from the literature, we show that, at fixed circular velocity, the stellar 

mass of star-forming galaxies has increased by a factor 2.5 between z = 2 and z — 0, 
whilst the rest-frame B-band luminosity has decreased by a factor ^ 6 over the same 
period. Together, these demonstrate a change in mass-to-light ratio in the i?-band 
of A(M /Lb) / (M/Lb)z=o "^3.5 between 2 = 1.5 and z — 0, with most of the evolu- 
^ ■ tion occuring below 2 = 1. We also use the spatial variation of [Nil] / Ha to show that 

I the metallicity of the ionised gas in these galaxies declines monotonically with galac- 

tocentric radius, with an average A log(0 / H) / A R= -0.027± 0.005 dexkpc"^ This 
gradient is consistent with predictions for high-redshift disk galaxies from cosmologi- 
cally based hydrodynamic simulations. 
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1 INTRODUCTION 

Numerical simulations indicate that the majority of mas- 
sive galaxies at z = 1-3, the epoch when galaxies were 
most rapidly growing their stellar mass, are continu- 
ously fe d by gas which promotes and maintains star for- 
mation ( jBournaud fc Elmcgrocnl |2009| : iDekel et al.1 l2009al : 
Ivan de Vo ort ct al. 2011, 2012 ). The accretion of gas from 
the halo, through minor mergers or by cold streams from the 
inter-galactic medium appears to dominate over that from 
major mergers and suggests that the star formation within 
the inter-stellar medium (ISM) of high-redshift galaxies are 



driven by internal dynamical processes iKeres et al.ll2005l ). 
The observational challenge is now to quantitatively mea- 
sure the internal properties (e.g. gas surface density, disk 
scaling relations, chemical abundances, distribution and in- 
tensity of star-formation) and so test whether the prescrip- 
tions developed to describe star-formation processes within 
disks at 2 = are applicable in t he rapidly evolving ISM of 
gas rich galaxies at h igh-redshift (|Krumholz fc McKeell2005l : 
iHopkins et aLll2012l 'l. 

High-spatial resolution observations of star-forming 
galaxies around z 1-2 have shown that a large fraction 
of the population have their ionised gas in large, rotating 
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disks (e.g. iGen zel et al. ll2006l: iForster Schreiber et al.ll2006l. 



20091: IStark et al. 200&r ijones et al.ll2010bl : IWisnioski et all 



20111 ). In contrast to low-redshift, these disks have high gas 



fractions (/gas = 20-80%) and so are turbulent, with much 
higher velocity dispersions given their rotational velocity 
(cr = 30-100 km s~^, Vrot / cr ~ 0.2-1) compared to the thin 



disks of loca,! spirals (iTacconi et al 



.20101 : iDaddi et al.ll2010l : 

iGeach et al]|201ll : [s^nbank et al.ll2011^ . 

If low- and high- redshift disks can be linked in a single 
evolutionary model, then the redshift evolution of disk scal- 
ing relations (i.e. between size, luminosity, rotational speed 
and gas fraction) providing a diagnostic of the relationship 
between the recent- and past averaged- star-formation with 

the dark halo and the domina nt mode of assembly. 

The pioneering work of IVoet et all l| 19961 ) and IVogtl 

l|l999l ) provided the first systematic study of the evolu- 
tion of the stellar luminosity (Mb) ve rsus circular veloc- 
ity (Tully-Fisher; ITuUv fc Fished 1 19771 ) relation, deriving 
an increase in luminosity at a fixed circular velocity of 
ALb ~ 1.7 X up to z ~ 1, reflecting the cha nge in star for- 
mation efficiency over this peri od ( see also Bamford et al. 
l2005l : lFernandez Lorenzo ct al. 2003 andlMiller et al. 2011). 
Since the B-band luminosity is sensitive to recent star 
formation, attempts have also been made to measure the 
evolution of the stellar mass (M*) Tully-Fisher relation 
which reflects the relation between the past-average star- 
formation history and halo mass. Although the current sam- 
ples sizes are drawn from a heterogeneous mix of popu- 
lations, evidence is accumulating that there is only mod- 
est evolution in the stellar-mass Tully-Fisher relation to 
z ~ 1 .5, with stronger evolution above z ~1.5 (jMiller et al.l 
I2OI2I ). Theoretical mode of galaxy formation have strug- 
gled to reproduce the zero-point and slope of the Tully- 
Fisher relation, tending to pr oduce dis ks that rotate too fast 
given their luminosities (e.g. IMo et a l. 1998; Bens on et al.l 
l2003l : iDutton et al]l201ll '). although some of this discrep- 
ancy has been alleviated by improving the recipes for star- 
burst driven feedback (e.g. Sales et al. 2010l: Tbnini et akl 



2011: 'Piontek & Steinmetz' 2011': 'Sc annapieco et al.l I2OIII : 
Desmond 2012; McCa rthy et al. 2012). 

A second potential observational consequence of differ- 
ences in fueling mechanisms for disk galaxies is their chemi- 
cal abundance distributions. The primary indicator of chem- 
ical evolution is typically traced by Oxygen as its relative 
abundance surpasses all elements heavier tha n Helium and 
it is present almost entirely in the gas phase (jSnow fc Witd 
ll996l : ISavage fc Sembachll 19961) . 

In the thin disk of the Milky- Way and other lo- 
cal spirals, t here is a ne g ative radial metallicity gra- 
dient (e.g 



Searlj Il97ll: IShieldj Il974l: iMcCall et all 



19851: IVila-Costas fc Edmundd 119921: IZaritskv et all Il994l: 



Ferguson et al ] |l998l :l van Zee et al.l[l998l ). In contrast, the 



thick dis k of the Milky Way displays no radial abundance 
gradient dCilmore et al.lll995l : [Beiilll996l : iRobin et all 1 19961 : 
lEdvardsson et all 1 19931 ). Thick disks are ubiquitous in spi- 
ral galaxies, and typically contain 10-25% of their baryonic 
mass, yet the thick di sk of the Milky- W ay contains no stars 
younger than 8 Gyr (jCilmore fc Wvsj Il9 85: Re ddv et al.l 
I2OO6I ). Thick disks are comparable in many physical proper- 
ties to the apparently turbulent disks of high-redshift galax- 
ies, and so determining the early form of their abundance 
gradient (positive or negative) will provide a key indica- 



tion of how quickly enrichment proceeds in inner and outer 
regions, as well as the interplay between star formation, 
clump migration and gas accretion from the halo or IGM 
( Solway et al. 2012). 

Indeed, if the majority of the gas accretion in high- 
redshift galaxies is via accretion from the IGM along fll- 
aments which intersect and deposit pristine material onto 
the galaxy disk at radii of 10-20 kpc, then the inner disks 
of galaxies should be enriched by star-formation and su- 
pernovae whilst the the outer-disk continually diluted by 
pristine gas, leavin g strong negative abundance gradients 
JPekel et al.ll2009al ) which will flatten as the gas accretion 
from the IGM becomes less efficient, and the gas redis- 
tributed, at lower redshift. 

To gain a complete census of how galaxies assemble the 
bulk of their stellar mass, it is clearly important to study the 
dynamics, chemical properties and internal star-formation 
processes in a well-selected sample of high-redshift galax- 
ies. To this end, we have combined panoramic (degree-scale) 
imaging with near-infrared narrow-band filters, to carry out 
the High-Z Emission Line Survey (HiZELS) survey, target- 
ing Ha emitting galaxies in four preci se (5z = 0.03) red- 
shift slices: z = 0.40, 0.84 1.47 fc 2.2 3 ijGeach et al.ll2008l : 
ISobral et al.l |2009| . |2010| . l201ll . l2012al lb[). This survey pro- 
vides a large, luminosity-limited sample of identically se- 
lected Hq emitters at epochs spanning the apparent peak of 
the cosmic star-formation rate density, and provides a pow- 
erful resource for studying the properties of representative 
galaxies (i.e. star-formation rates ~lO'sM0/yr) which will 
likely evolve into ~ L* galaxies by z — 0, but are seen at a 
time when they are assembling the bulk of their stellar mass, 
and thus at a critical stage in their evolutionary history. One 
of the key benefits of selecting a high-redshift galaxy pop- 
ulation using Hq alone is the simplicity in comparison to 
similar luminosity-limited samples of star-forming galaxies 
in the local Universe. Together, these comparison samples 
can be used to search for differences in the physical prop- 
erties of the ISM that potentially provide insight into the 
factors infiuencing star formation at high redshift. 

In this paper, we present adaptive optical assisted 
(~kpc scale) integral field spectroscopy with SINFONI of 
nine star-forming galaxies selected from the HiZELS survey 
in the redshift range z = 0.84-2.23 - the SINFONLHiZELS 
(SHiZELS) survey. We use these data to investigate the 
dynamical properties of the galaxies, the evolution of the 
luminosity and stellar mass scaling relations (through the 
Tully-Fisher relation), and the star formation and enrich- 
ment within their ISM. We use a cosmology with Qa ~ 0.73, 
flm = 0-27, and Ho = 72 kms~^ Mpc~^. In this cosmology, at 
the median redshift of our survey, 2=1.47, a spatial res- 
olution of 0.1" corresponds to a physical scale of 0.8 kpc. 
All quoted ma gnitudes are on the AB system and we use a 
Chabrier IMF (jchabrieill2003l ). 



2 SAMPLE SELECTION, OBSERVATIONS & 
DATA REDUCTION 

2.1 HiZELS 

To select the targets for IFU observations, we exploited the 
high density of sources in the HiZELS imaging of the COS- 
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Figure 1. Broad-band photometry and SEDs for the nine galaxies in our sample from rest- frame- UV to mid-infrared bands (spanning 
GALEX far-UV and near-UV bands to Spitzer /IRAC) . In all cases, we calculate the rest-frame spectral energy distribution (SED) and 
sta r-formation history and use this information to estimate the stellar masses (Table 2). The SEDs are calculated using a fit using 
the lBruzual fc Chariot {2003) population synthesis models with a Chabrier IMF. The SHiZELS galaxies sample a range of stellar masses 
from ~10^~^^ Mq and specific star-formation rates 0.1— lOGyr"^. 



Table 1: Targets L Observations 



ID 




RA 


Dec 


^Hn 


tcxp 


Strchl 


EE 






(;J2000) 


(.12000) 






(ks) 




(0.1") 


SHiZELS 


1 


02 18 26.3 


-044701.6 





8425 


9.8 


15% 


11% 


SHiZELS 


4 


10 01 55.3 


+ 02 14 02.6 





8317 


3.6 


10% 


S% 


SHiZELS 


7 


02 1700.4 


-05 01 50.8 


1 


4550 


13.4 


27% 


25% 


SHiZELS 


S 


02 IS 21.0 


-05 19 07.8 


1 


4608 


9.8 


21% 


27% 


SHiZELS 


9 


02 17 13.0 


-04 54 40.7 


1 


4625 


9.8 


24% 


27% 


SHiZELS 


10 


02 1739.0 


-04 44 43.1 


1 


4471 


9.8 


20% 


25% 


SHiZELS 


11 


02 IS 21.2 


-05 02 4S.9 


1 


4858 


3.6 


13% 


32% 


SHiZELS 


12 


02 19 01.4 


-04 58 14.6 


1 


4676 


9.8 


20% 


30% 


SHiZELS 


14 


10 00 51.6 


+02:33 34.5 


2 


2418 


12.0 


12% 


26% 



Table 1. Notes: EE denotes the encircled energy within a radius 
of 0.1". 



MOS and UDS fields l|Sobral et al.ll2012bl ') to select Ha emit- 
ters which are close to bright (R < 15.0) stars, such that nat- 
ural guide star adaptive optics correction could be applied 
to achieve high spatial resolution. For this program, we se- 



lected fourteen galaxies with Ha fluxes from the HiZELS 
catalog with fluxes between 0.7-1.6 x 10"^^ ergs~^ cm~^ to 
ensure that their star-formation properties and dynamics 
could be mapped in a few hours. However, of the fourteen 
galaxies in our parent sample, only nine systems yield a Ha 
detection in the data-cubes, which may (in part) be due to 
spurious detections in the early HiZELS catalogs. At the 
three redshift ranges of our sample, the average Ha fluxes 
of our galax ies correspond to star-formation rates (adopt- 
ing the the iKennicutO 1 19981 calibration with a Chabrier 
IMF) of SFR = 4.6 x lO"** Lhq = 1.7, 6.9 and 19.3 Mq yr'^ 
at 2 = 0.83, 1.47 and 2.23 respectively. 

Next we estimate the far-infrared luminosities for the 
galaxies in our sample. We exploit the Herschel/SPTRE 
imaging of these flelds which were taken as part of the 
Herschel Multi-tiered Extra-galactic Survey (HerMES). The 
SPIRE maps reach 3 a limits of ~8, 10 and llmjy at 250, 
350 and 500/im respectively. None of the SINFONI-HiZELS 
(SHiZELS) galaxies in our sample are individually detected 
(which is unsurprising given the implied star-formation rates 
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Table 2: Integrated Galaxy Properties 



ID 


ergs"-*- cm~^) 


SFRhq 
(Mo/yr) 


'"1/2 

(kpc) 


[Nil] / Ho 


Ms 
(AB) 


Mk 
(AB) 


logio(M.) 
(M©) 


E (B - V) 


A log(0 / H) / A R) 
(dex kpc~^ ) 


SHiZELS-1 


l.OitO.l 


2 


1.8±0.3 


0.2 ±0.05 


-18.70 


-21.54 


10.03 ±0.15 


0.4±0.1 


-0 037+0-''^'' 


SHiZELS-4 


0.7±0.1 


1 


1.4±0.5 


0.2±0.1 


-20.68 


-21.49 


9.74 ±0.12 


0.0±0.2 


SHiZELS-7 


1.2±0.1 


8 


3.7±0.2 


0.43 ±0.05 


-21.12 


-22.15 


9.81 ±0.28 


0.2 ±0.2 


-0 019+°-"^^ 


SHiZELS-8 


l.litO.l 


7 


a.iito.s 


<0.1 


-21.26 


-21.90 


10.32 ±0.28 


0.2 ±0.2 


+0 006+° ''^^ 


SHiZELS-9 


l.OitO.l 


6 


4.1±0.2 


0.27±0.03 


-21.77 


-22.49 


10.08 ±0.28 


0.2 ±0.2 


-0 027+°"^" 
'J-'^^' -0.018 


SHiZELS-10 


1.6±0.2 


10 


2.3±0.2 


0.13 ±0.04 


-21.15 


-21.81 


9.42 ±0.33 


0.3 ±0.2 


-0 031+°"^'^ 


SHiZELS-11 


1.3±0.1 


8 


1.3±0.4 


0.6±0.1 


-22.13 


-24.44 


11.01 ±0.24 


0.5 ±0.2 


-0.0870tO;«g2 


SHiZELS-12 


0.7±0.1 


5 


0.9±0.5 


0.27±0.03 


-21.34 


-22.44 


10.59 ±0.30 


0.3±0.2 


SHiZELS-14 


1.6±0.1 


27 


4.6±0.4 


0.60 ±0.05 


-23.55 


-24.75 


10.90 ±0.20 


0.4±0.1 


-0 024+° ''^2 


Median 


l.litO.l 


7±2 


2.4 ±0.7 


0.27 ±0.08 


-21.3±0.20 


-22.2±0.3 


10.25 ±0.50 


0.3±0.1 


-0.026 ±0.006 



Table 2. Notes: Ha fluxes are in units of 10~^® ergs~^ cm~^ . The errors on the fluxes are dominated by ~10% systematic uncertainties 
in the flux calibration from the change in Strehl ratio between observing standard stars on-axis and the tar get galaxies o ff axis. r]^y2 is 
the Ha half light radius and has been deconvolved for the PSF. Errors on Mb & Mk are typically 0.15 dex llSobral et al.ll201ll) . 



from Ha), but we can search for an average star- formation 
rate via stacking. 

In the following analysis, we restrict ourselves to the 
z = 1.47 galaxies only to mitigate against redshift effects in 
the stack on the SED. The 250, 350 & 500^111 stacked fluxes 
for the 2 = 1.47 SHiZELS galaxies of 6.9 ±1.7, 9.0 ±1.9 
and 7.0±2.4mJy at 250, 350 and 500/im respectively (or 
7.5 ±1.6, 9.2 ±2.2 & 7.4±3.0mJy at 250, 350 & 500/.im if 
we instead stack all nine SHiZELS galaxies in our sample). 
A simple modified black-body fit to this photometry (with a 
dust emmisivity of /3 = 1.5-2.0) at the known redshift yields 
a far-infrared luminosity (integrated between 8-1000/.im) of 
Lf/h = 1.8 ± 0.6 X 10^^ L© and hence a total star-formation 
rate of SFRfir = 2.7 x 10"'"' Lfir = 18 ± 8 M© yr'V 

We note that if we instead stack all 317 HiZELS 
Ha emitters from the z=1.47 parent sample with 
fluxes >1 X 10"^®ergs"^cm~2 in the UDS and COS- 
MOS fields then the 250, 350 & 500/.im fluxes are 
6.3 ±0.24, 7.9 ±0.25 and 5.6 ± 0.31 mjy respectively, sug- 
gesting LpiR = 1.5 X 10^^ L© which is consistent with our 
SHiZELS sample. 

Next, for all of the SHiZELS galaxies, we use the 
broad-band imaging of the galaxies in order to calculate 
the rest-frame spectral energy distribution (SE D), reddening 
and s tar-formation histories and stellar masses (|Sobral et al.l 
l20ld ). We exploit deep archival multi-wavele ngth imaging 
from rest-frame UV to mid-infrare d bands (see Sobral et al.l 
I2OIOI for references) and follow ISobral et aO (120111) to 
perfor m a full SED fit using the 'Bruz ual fc CharlotI 
l|2003l ) and Bruzual (2007) population synthesis models 
which include the Thermally Pulsating Asymptotic Giant 
Branch (TP-AGB ) stellar phase and a Chabrier IMF (see 
ISobral et al]|2010l for the range of parameters) . We use pho- 
tometry from up to 36 (COSMOS) and 16 (UDS) wide, 
medium and narrow bands (spanning GALEX far-UV and 
near-UV bands to Spiteer/IRAC; Fig[l| and in Table 2 we 
give the best-fit stellar mass and E(B— V) for each galaxy. 

The average E(B— V) for our sample is 
E(B-V) = 0.3 ± 0.1 which corresponds to Ahc = 0.91 ± 0.21 



(Av = 1.11 ± 0.27). Applying this correction to the Ha 
luminosities suggests reddening corrected star-formation 
rates SFRhq = 16 ± 5 M© yr~^, which is similar to those 
derived from the far-infrared (SFRfir = 18 ± 8 M© yr~^). 

The rest-frame B-band magnitudes (uncorrected for ex- 
tinction) and stellar masses are given in Table 2 (the Icr 
errors are drawn from the multidimensional distribution 
and are typically 0.15 dex). So that more direct comparisons 
can be made with other datasets, we also report the rest- 
frame B- and A'-band magnitudes for each of the galaxies 
in our sample, and we note that the average inferred light- 
to-mass ratio is Lk/M = 0.58 ± 0.14. 

Thus, overall, these galaxies appear to have luminosi- 
ties consistent with local LIRGs, with dust corrections of 
Av ^ 1 and star-formation rates char acteristic of "typical " 
star- forming galaxies at these epochs jSobral et al.|[2012bh . 
Moreover, the SHiZELS galaxies sample a range of stellar 
masses from ~10'~^^ M© and specific star-formation rates 
0.1-10 Gyr~^. Since one of the main high-redshift compari- 
son samples we use throughout thi s paper is the SINFONI 
Near -Infrared Galaxy Survey (SINS: lForster Schreiber et al.l 
i200^, we note that our sample has a comparable range of 
stellar masses, although the SINS sample tend to have higher 
Ha-derived star-formation rates, (SFRhu = 7 ± 2 M© yr~^ 
and 3O±3M0yr~^ for SHiZELS and SINS respectively). 

2.2 SINFONI Observations 

To map the nebular Ha emission line properties of the galax- 
ies in our sample, we used th e SINFONI IFU on the ESO 
VLT (|Eisenhauer^ et al.ll2003l ). The SINFONI IFU uses an 
image slicer and mirrors to reformat a field of 3 x 3" at a 
spatial resolution of 0.l"/pixel. At 2 = 0.84, 1.47 and 2.23 
the Ha emission line is redshifted to ~ 1.21, 1.61 and 2.12/im 
and into the J, H and A-bands respectively. The spectral 
resolution in each band is A / AA ~ 4500 (the sky lines have 
4- A FWHM). In the remainder of our analysis, emission 
line widths are deconvolved for the instrumental resolution. 
Since each of the targets were selected to be close to a bright 
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AO star, we used NGS correction. The -R-band magnitude 
of the AO stars are R = 12.0-15.0 and they have separations 
between 8" and 40" from the target galaxies. 

To observe the targets we used ABBA chop sequences, 
nodding 1.5" across the IFU. We observed each target for 
between 3.6 and 13.4 ks (each individual exposure was 600 
seconds) between 2009 September 10 and 2011 April 30 in 
~0.6" seeing and photometric conditions. The median Strehl 
achieved for our observations is 20% (Table 1) and the me- 
dian encircled energy within 0.1" is 25% (the approximate 
spatial resolution is 0.1" FWHM or 850 pc at 2 = 1.47 - the 
median redshift of our survey) . Individual exposures were re- 
duced using the SINFONl esorex data reduction pipeline 
which extracts, flat-fields, wavelength calibrates and forms 
the data-cube for each exposure. The final data-cube was 
generated by aligning the individual data-cubes and then 
combined these using an average with a 3-o" clip to reject 
cosmic rays. For flux calibration, standard stars were ob- 
served each night either immediately before or after the sci- 
ence exposures. These were reduced in an identical manner 
to the science observations. 



2.3 Galaxy Dynamics 



All nine galaxies display strong Ha emission, with a range of 
Ha luminosities of Lhci ~ lO*^'*"*^ * ergs~^. In Figure[2]we 
show the galaxy-integrated one-dimensional spectrum from 
each galaxy. In all cases we clearly detect the [Nn]A6583 
emission, the median [Nil] / Hq for the sample is 0.3 ±0.1, 
with a range of 0.1 < [Nil] / Ha< 0.6. None of the galaxies 
display strong AGN signatures in their near-infrared spec- 
tra. 

To search for fainter lines we also transform each of 
the spectra to the rest-frame and coadd them (weighted 
by flux) and show this stack in Fig. [2] In this compos- 
ite, the [Nil] /Ha = 0.29 ±0.02, is consistent with the av- 
erage of individual measurements. We also make a clear de- 
tection of the [Sii]AA 6716,6761 doublet and derive a flux 
ratio of lerie / lerai = 1-36 ± 0.08 which correspo nds to an 
electr on density in the range 100-1000 cm~^ fjOsterbrockl 
1 19891 ) and a mass in the ISM of 3-30 x 10^" Mq for a half- 
light radius of 2kpc (Table 2). It is also useful to note 
that the [Sll] A6716 / Ha reflects the ionisation strength of 
the ISM and we derive [Sll] A6716 / Ha = 0.13 ± 0.04 which 
suggests an ionisation parameter log(U /cm^) = —3.5 ± 0.3 
(|Osterbrocklll989l : ICollins fc Randll200ll ). This is compara- 
ble to an ionisation rate of logjQ(U / cm^) ~ —3.9 which can 
be inferred from the average star-formation rate volume den- 
sity of our galaxies and assuming there ar e 3 x 10^° ionising 
photons per solar mass of star-formation (|Coxlll983l '). 

To measure the velocity structure in each galaxy, we fit 
the Ha and [Nil] AA6548, 6583 emission lines pixel- by-pixel. 
We use a minimisation procedure, taking into account the 
greater noise at the positions of the sky lines. We first at- 
tempt to identify a line in each 0.1 x 0.1" pixel (~ 1 x 1 kpc, 
which corresponds to the approximate PSF), and if the fit 
fails to detect the emission line, the region is increased to 
2x2 pixels. Using a continuum fit we required a signal-to- 
noise >5 to detect the emission line in each pixel, and when 
this criterion is met we fit the Ha and [Nil] AA 6548,6583 
emission allowing the centroid, intensity and width of the 
Gaussian profile to vary (the FWHM of the Ha and [Nil] 



lines are coupled in the fit). Uncertainties on each param- 
eter of the fit are calculated by perturbing one parameter 
at a time, allowing the remaining parameters to find their 
optimum values, until Ax^ = 1 is reached. 

Even with ~kpc-scale resolution, we note that there 
is a contribution to the line widths of each pixel from the 

large-scale velocity motions across the PSF which must be 

I ^ II — ' — I 

corrected for (jPavie s ct al. 2011'). To account for this "beam 

smearing" , at each pixel where the Ha emission is detected, 
we calculate the local velocity gradient (A V A R) and sub- 
tract this in quadrature from the measured velocity disper- 
sion. 

In Fig. [3] we show the Ha intensity, velocity and line of 
sight velocity dispersion maps for the galaxies in our sample. 
There is clearly a variety of Ha morphologies, ranging from 
compact (e.g. SHiZELSll & 12) to very extended/clumpy 
(e.g. SHiZELS 7, 8, 9, 14). For three galaxies, we also include 
high-resolution broad-band imaging thumbnails from Hub- 
ble Space Telescope ACS V-hand (SHiZELS 4&14) or WFC3 
_ff-band (SHiZELS 11) which shows that the Ha is well cor- 
related with the rest-frame UV/optical morphology. The dis- 
tribution and intensity of star-formation (and the properties 
of the star-forming clumps) are discussed in a companion pa- 
per (Swinbank et al. 2012 ApJ submitted). Fig.[3]also shows 
that there are strong velocity gradients in many cases (e.g. 
SHiZELS 1, 7, 8, 9, 10, 11, & 12). We will investigate how 
many of these systems can be classified as "disks" in §3. 



3 ANALYSIS, RESULTS & DISCUSSION 

All nine galaxies display velocity gradients in their Ha kine- 
matics, with peak-to-peak velocity gradients ranging from 
40-400 kms~^ and a ratio of peak-to-peak difference (umax), 
to line of sight velocity dispersion (a) of -Umax / o' — 0.3-2.9. 
Although a ratio of «max / o- — 0.4 has traditionally been 
use d to crudely differentiate rot ating systems from merg- 
ers (|Forster Schreiber et al.ll2009l ). more detailed kinematic 
modelling is essential to reliably differentiate rotating sys- 
tems from mergers (jSliapiroj5t^al. , 200^). We therefore begin 
by modelling each velocity field with a rotating disk in order 
to estimate the disk inclination and true rotational velocity 
(we exclude SHiZELS 4 & 12 from this analysis as their ve- 
locity fields are only marginally resolved by our data). 

Although there have been several functional forms 
to descri be rotation curv es (such as the "multi-parameter 
Courteaul Il997l or the "universal rotation curve; 



fit": e.g.. 

iPersic et al.l Il996l ). the simplest model which provides a 
good fit to most rotation curves on ~ kpc scales is the arc- 



tan fit of the form v{r) -- 



1 arctan(r/ri), where «asym is 



the asymptotic rotational velocity and rt is the effective ra - 
dius at which the rotation curve turns over (Courteaulll997l). 
To fit the two-dimensional velocity fields of the galaxies in 
our sample, we construct a two-dimensional kinematic model 
with six free parameters (vasym, rt, [x/y] centre, position an- 
gle (PA) and disk in clination) and use a genetic algorithm 
(|Charbonne"aulll995D to find the best model. During the fit- 
ting, we limit the parameters ranges to those allowed by the 
data (e.g. the [x/y] center and turn over radius, rt, must 
be within the SINFONl field of view; the maximum cir- 
cular velocity must be < 2 x the maximum velocity seen 
in the data; PA = 0-180°, and i = 0-90). We then generate 
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Figure 2. Spatially integrated, one dimensional spectra around the redshifted Ho emission for each of the galaxies in our sample. In 
each case, the dashed lines identify the Ha, [Nil] AA 6583,6548 and [Sll]AA 6716,6731 emission lines. In the final panel we show the flux 
weighted composite spectrum. In all cases, we clearly detect [Nil] emission and the median [Nil] /H« for the sample is 0.3 ±0.1, with 
a range of 0.1 < [Nil] / Ha < 0.6. None of the galaxies display strong AGN signatures in their near-infrared spectra (e.g. broad lines or 
strong [Nil] / Ha ratios). 



Table 3: Dynamical Properties 



ID 


Vasym 


C7 


V2.2 


Inc 


Kv 




KTot 




Class 




(km/s) 


(km/s) 


(km/s) 


(deg) 












SHiZELS-1 


112±11 


98±15 


86±18 


40 


0.12±0.17 


0.17±0.81 


0.20±1.03 


4.8 


D 


SHiZELS-4 




77±20 














C 


SHiZELS-7 


145±10 


75±11 


112±15 


33 


0.10±0.08 


0.35±0.18 


0.36±0.34 


3.1 


D 


SHiZELS-8 


160±12 


69±10 


144±17 


25 


0.20±0.15 


0.30±0.25 


0.36±0.40 


3.0 


D 


SHiZELS-9 


190±20 


62±11 


155±24 


54 


0.12±0.05 


0.48±0.20 


0.49±0.30 


5.7 


D 


SHiZELS-10 


30±12 


64±8 


26±12 


30 


0.13±0.05 


0.83±0.14 


0.84±0.35 


25.5 


M 


SHiZELS-11 


224±15 


190±18 


208±18 


25 


0.16±0.05 


0.34±0.10 


0.37±0.16 


5.34 


D 


SHiZELS-12 




115±10 














C 


SHiZELS-14 




131±17 


190±25 


58 


0.15±0.09 


0.38±0.23 


0.40±0.34 


9.1 


M 


Median 


147±31 


75±19 


133±30 




0.13±0.01 


0.37±0.07 


0.38±0.04 


5.3±3.0 





Table 3. Notes: Vasym and V2.2 are inclination corrected, cr is the avergae velocity dispersion across the galaxy image (and the velocity 
dispersion map has been corrected for "beam smearing" effects due to the PSF, see § refsec:dynamics. The classes in the final column 
denote Disk (D), Merger (M) and Compact (C) (see §3) 



10* models within random initial parameters between these 
ranges and find the lowest xi- This process is repeated (up 
to 100 times), but in each new generation we remove (up to) 
10% of the worst fits and contract the remaining parameter 
space search by up to 5% centered on the best-fit. We ex- 
amine whether the fit has converged by examining whether 
all of the models in a generation are within Ax^ = 1 of the 
best-fit model. The errors on the final parameters (Table 3) 
reflect the range of acceptable models from all of the models 
attempted. 

The best-fit kinematic maps and velocity residu- 
als are shown in Fig. O whilst the best-fit inclina- 
tions and disk rotation speeds are given in Table 2, 
along with the x2- AH galaxies show small-scale devia- 
tions from the best-fit model, as indicated by the typi- 
cal r.m.s, < data — model > = 30 ± 10 kms~^, with a range 
from < data — model > = 15-70 km s^^. These offsets could 
be caused by the effects of gravitational instability, or sim- 
ply due to the unrelexed dynamical state indicated by the 
high velocity dispersions (ct = 96 ± 24 kms~^). The goodness 



of fit and small scale deviations from the best-fit models are 
similar to those of 18 galaxies in the SINS survey which 
show the most prominent rotational support (r.m.s of 10- 
80kms"^ and =0.2-20; ICresci et aLl(2009l ). as well as to 
the hig h-resolution studies of gr avitationally lense d galax- 
ies from I Jone s et al.l l|2010bl ) and lStark et alj (|2008f) (which 
have an average x^ =3.5 ±0.8). We therefore conclude that 
the disk model provides an adequate fit to the majority of 
these galaxies and that the velocity fields of most galax- 
ies are consistent with the kinematics of turbulent rotating 
disks. 

We use the best-fit dynamical model to identify the ma- 
jor kinematic axis, and extract the one dimensional rota- 
tion curve and velocity dispersion profiles, and show these 
in Fig. 131 We define the zero point in the velocity using the 
dynamical center of the galaxy, whilst the error bars for the 
velocities are derived from the formal Icr uncertainty in the 
velocity arising from the Gaussian profile fits to the Ha emis- 
sion. We note that in the one-dimensional rotation curves in 
Fig. [4] alternate points show independent data. We also note 
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Figure 3. Ha and kinematic maps of the SHiZELS galaxies. For each galaxy, the left hand panel shows the Ha emission line flux. The 
contours denote a Ha-derived star-formation surface density of Ssp = 0.1 Mq yr~^ kpc"'^ . The central two panels show the velocity field 
and line-of-sight velocity dispersion profile (tr) respectively. The right hand panel shows the residual velocity field after subtracting the 
best-fit kinematic model. The r.m.s. of these residuals is given in each panel (for SHiZELS 4 & 12 there are too few resolution elements 
across the source to meaningfully attempt to fit disk models). 



that the data have not been folded about the zero velocity, 
so that the degree of symmetry can be assessed. 

Whilst this modelling provides a useful test of whether 
the dynamics can be described by disk-like rotation, a better 
criterion for distinguishing between motion from disturbed 
kinematics is the "kinemetry" which measures the asymme- 
try of the velocity field and line-of-sight velocity dispersion 
maps for each galaxy. This technique has been well cal- 



ibrated and tested at low redshift (e.g. iKrainovic et al.l 
[2006r), whilst at high redshift it has been used to deter- 
mine the strength of deviations of the observe d velocity and 
dispe rsion maps from an ideal rotating disk l|Shapiro et al.l 
|2008| ). In this modelling, the velocity and velocity disper- 
sion maps are described by a series of concentric ellipses of 
increasing semi-major axis length, as defined by the system 
center, position angle and inclination. Along each ellipse, 
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Figure 3. continued. 



the moment map as a function of angle is extracted and de- 
composed into it s Fourier series which have coefficients fcn 
at each radii fsee lKrainovic et al.ll2006l for more details). 



In a noiseless, idealised disk, the coefficients of the 
Fourier expansion would be zero. Any deviations from the 
ideal case that occur in a disturbed disk (e.g. warps, spi- 
ral structures, streaming motions or mergers in the extreme 
cases) induce strong variations in the higher-order Fourier 
coefficients. To ensure that reliable measurements of the 
kinemetry parameters can be made, we first construct a grid 
of 2000 model disks and mergers with noise and spatial res- 
olution appropriate for our SINFONI observations. For the 
model disks, we allow disk scale lengths and rotation curve 
turn over lengths (rt) from 0.5-3 kpc, rotation speeds and 
peak velocity dispersions from 10-300 km s"'^ and random 
inclinations and sky position angles. The model mergers typ- 
ically comprise 2-3 of these galaxy-scale components sepa- 
rated by 1-20 kpc in projection. The dynamics of each of 



these components randomly ranges from dispersion domi- 
nated to rotating systems. 

W e then use the kinemetry code from iKrainovic et al.l 
(|2006h and measure the velocity and veloc ity dispersion 
asymm etry for all of these models. Following IShapiro et al.l 
(|2008l ). we define the velocity asymmetry (Kv) as the av- 
erage of the kn coefficients with n = 2-5, normalised to the 
first Cosine term in the Fourier series (which represents cir- 
cular motion); and the velocity dispersion asymmetry (Kcr) 
as the average of the first five coefficients (n = 1-5) also nor- 
malised to the first Cosine term. For an ideal disk, K„ and 
Ka will be zero. In a merging system, strong deviations from 
the idealised case causes large Ki, and values, which can 
reach Kv ~ Ko- ~ 10 for very disturbed systems. The total 
asymmetry, Ktoi is -ft^xot = + . For our mock sample 
of model disks, we recover Kxot.disk = 0.30 ± 0.03 compared 
to Kxot, merger = 2 ± 1 fot the mcrgcrs. 

For the SHiZELS galaxies in our sample, we measure 
the velocity and velocity dispersion asymmetry, (SHiZELS 4 
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Figure 4. Left: One dimensional rotation curves for tlie seven 
galaxies in our sample which show strong velocity gradients in 
their Ha dynamics. For the six galaxies whose dynamics resemble 
rotating systems, we overlay the best-fit one dimensional rotation 
curve Right: One dimensional velocity dispersion profiles with the 
best-fit line of sight profile overlaid for the rotating systems. The 
observed rotation curve and line-of-sight velocity dispersion are 
extracted from the two-dimensional velocity fields using a 1 kpc 
width slit across the major kinematic axis of each galaxy. 



& 12 have too few independent spatial resolution elements 
acros s the galaxy so w e omit these from the kinemetry anal- 
ysis). iKrainovic et al ] (l20Q6l) show that an incorrect choice 
of centre induces artificial power in the derived kinemetry 
coefficients. We therefore allow the dynamical center to vary 
over the range allowed by the family of best-fit two dimen- 
sional models and measure the kinemetry in each case. We 
also perturb the velocity and dispersion maps by the errors 
on each pixel and re-measure the asymmetry, reporting the 
velocity and dispersion asymmetries, (Ky and Ko- respec- 
tively) along with their errors in Table 3. We use the total 
asymmetry, Ktoi to differentiate disks (Kiot < 0.5) from 
mergers Ktoi >0.5). For the galaxies in our sample, five 
have are classified as disks (D), whilst two more indicate 
mergers (M), and the final two are compact (C). 

As expected, the two galaxies that are classified as 
mergers from the kinemetry (SHiZELS 10 & 14) have the 
highest iXv > 9) from the dynamical modelling. Both of 
these systems are noteworthy; the one-dimensional rotation 
curve of SHiZELS 10 (extracted from the major kinematic 
axis) appears regular and suggests rotation, but the two- 
dimensional velocity field reveals extended (~ 5 kpc) emis- 
sion to the east, which may be due to a companion or tidal 
material. The dynamics of the other galaxy classified as a 
merger, SHiZELS 14, are clearly more complex. The Ha is 
extended across ~ 14kpc with a peak-to-peak velocity gra- 
dient of ~480 ± 40 kms~^. Qualitatively, the one- and two- 
dimensional velocity field of SHiZELS 14 are consistent with 
an early stage prograde encounter, with each component in 
the merger displaying approximately equal peak-to-peak ve- 
locity gradients. The northern most component is also more 
highly obscured than the southern component (as evident 
from the HST 7-band imaging; Fig. |3J . 

Overall, the fraction of rotating systems within our 
sample, ~ 55%, is consistent with that found from other 
Ha IF U surveys of high-red s hift st ar-forming galaxies 



Forster Schreiber et all l2009l : Ijones et al.l l2010bl : 



IWisnioski et al.ll201ll ) 



3.1 The Tully-Fisher Relation 

Since our two dimensional dynamical data resolve the 
turn over in the rotation curves (as evident from Fig. 3}, 
we can use our results to investigate how the disk scal- 
ing relations for the galaxies in our sample compare to 
galaxy disks at z = 0. The relation between the rest- 
frame B-band luminosity (Mb) and rotational velocity and 
that between the total stellar-mass (M*) and rotational 
veloc ity define the Tully-Fisher relations ||Tu11v fc Fisher! 
[i973). The first of these relations has a strong contribu- 
tion from the short-term star-formation activity whilst the 
second provides a better proxy for the integrated star- 
formation history. Indeed the latter relationship may re- 
flect how rotationally-supported galaxies formed, perhaps 
suggesting the presence of self-regulating processes for star- 
formation in galactic disks. The slope, intercept and scatter 
of the Tully-Fisher relations and their evolution are there- 
fore key parameters that a ny successful g alaxy for mation 



model must reproduce (e . g. ICole fc Kaisei 



1989; 



|1994 iKauffmann et al] Il993l: lEisenstein fc Loebl 



Somerville fc Primackl 1 19991: [Ste inmctz fc Navarro! Il999l : 
IPortinari fc Sommer-Larsenll2007riDutton et aLlbOllI ). 
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Since we wish to compare our dynamical measurements 
with those of other high-redshift galaxies where similar data 
exist, it is important to measure the rotational velocity in a 
consistent manner. Whilst the arctan function is useful for 
describing the overall shape of the rotation curve (in both 
one and two-dimensional data), it provides an extrapolation 
of the rotational velocity, which may not be warranted es- 
pecially if the rotation curve does not flatten in the data. 
Equally, since any measurement of the observed peak-to- 
peak rotational velocity (Vmax) is sensitive to the radius at 
which it has been possible to detect the dynamics, this mea- 
surement varies from both object-to-object depending on 
both the extent of the galaxy and signal-to-noise of the ob- 
servations and can therefore give biased measurements when 
comparing samples. 

A more robust measurement is that of the velocity at 
2.2 times the disk scale length (V2.2). This has a better 
physical basis since it meas ures the peak rotational velocity 
for a pure exponential disk (|Freemanlll970l :[ Courteau fc 
Il997f ) and provides the clo sest match the r adio (21 cm) line 
widths for local galaxies (ICourteau|[l997l ). Of course, the 
light profile of high-redshift galaxies are unlikely to follow 
pure exponential disk profiles, but it provides a reasonable 
estimate of where the rotation curve flattens, and as we will 



see below, is a good approximation of both Vn 



and Va 



in our data. To measure V2.2 in our data, we therefore mea- 
sure the Ha intensity as a function of radius in our galaxies 
(centered at the dynamical center and using the inclination 
and position angle set by the best fit two-dimensional kine- 
matic model) and infer the disk scale length (which ranges 
from 0.7-2.7 kpc and has a median of 1.8 ± 0.7 kpc). We then 
measure the rotational velocity at 2.2 disk scale lengths and 
report these in Table 3. The errors on these quantities are 
derived from both the errors on the observed velocity at 
this radius and the formal errors on the two-dimensional 
model velocity fields. For reference, we note that owing to 
the high resolution and signal-to-noise of our data, we have 
been able to resolve the fiattening in the rotation curves, 
the ratio of the maximum (observed) peak-to-peak velocity 
(Vmax) to V2.2 in our sample is Vmax / V2.2 — 0.93 ± 0.05 
whilst the ratio of asymptotic velocity (Vasym) to V2.2 is 
VasH™/V2.2= 1.25 ±0.08. 

3.1.1 B-Band Tully-Fisher Relation 

In Fig. [5] we s how the z — Tully-Fisher relation from 
ITuUv fc Fisheil (|l977l ') and overlay the 2 = 0.84-2.23 mea- 
surements from our sample (using V2.2 as the preferred mea- 
sured of rotational veloci ty). We also inc l ude th e field sam- 
ple (< 2 > = 0.33) from iBamford et all (120051). the z: ~1 
field samples from iMiller et al.l (120111 ) a nd Iweiner et aT 
2006 ) , th e z = 1-3 lensed sa mples from ISwinbank et al 



20061) and lJones et al.l |2010bl') and the z ~ 2 ari d z ~3 star- 



forming galaxies from iForster Schreiber et al.l (|2009l ) and 
iGnerucci et al.l (|201ll ') respectively. 

It is clear from Fig.[5]that the z = _B-band Tully-Fisher 
relation provides a poor fit to high-redshift data, with most 
of the high-redshift galaxies lying above the local relation. 
Since the high-redshift sample represents a heterogeneous 
mix of populations, all with different and complex selection 
functions, we limit our search for the evolution in the Tully- 
Fisher to its zero-point only. We therefore fix the slope to 



the z — relation and in each of the seven redshift bins from 
z = 0.0-3 we derive the zero-point offset using a non-linear 
least squares regression (weighted by the velocity errors in 
each galaxy) and show these results in Fig. [5] 

For a fixed slope, the strongest evolution oc- 
curs between 2 = 0-1, and in this range the zero- 
point evolution follows AA/B(ma(7) = ( — 1.9 ± 0.3) 2 up to 
2 ~ 1 and then significantly flattens above 2 = 1 to 
AMb (ma5) = (-0.06 ± 0.25) 2-(l. 68 ±0.25) up to 2 ~ 3. 
Since the rest-frame _B-band is dominated by recent star- 
formation activity, the strong evolution in the rest-frame 
_B-band Tully-Fisher zero-point between 2 = and 2 = 1 
(2 magnitudes or a factor ~ 6 in luminosity) and then flat- 
tening to 2 ~3, is consistent with the e volution of the star- 
form ation volume density with redshift (|Hopkins fc BeacomI 
I2OO6I ). 

We can compare this evolution to galaxy formation 
models, and in Fig.[S]we also overlay the predic ted zero-point 
offset from the semi-an a lytic model s from Dutton et al.l 
(|201ll ) and iBower et al] (|2006l '). The iDutton et all (|201ll ) 
semi-analytic models suggest that for a fixed circular veloc- 
ity, the luminosity of a galaxy at 2 ~ 1 should be a factor up 
to ~2.5 X higher than today (|Dutton et al.|[201ll ). although 
the data are suggestive of markedly stronger evolution than 
predicte d by the models 2 ~ 1. In contrast the predictions 
from the lBower et al. I Hoo^) model predicts very little evo- 
lution in the i?-band Tully-Fisher relation {AMb < 0.2 mag- 
nitudes between 2 = 0-4). In this model, the apparent lack of 
evolution in the B-band Tully-Fisher relation in the galaxy 
formation models occurs because the specific star-formation 
rate evolves rapidly with both mass and redshift. Although 
the model approximately matches the evolution of the star- 
formation rate density with redshift, this is offset by an in- 
crease in the specific star-formation rate in low mass galax- 
ies. Moreover, since the dark halos are more concentrated 
at high-redshift, a fixed Vdrc does not select the same mass 
halo at the two epochs (for a fixed Vcirc, a halo at 2 = is 
a factor ~1.7x more massive than at 2 = 2). This combina- 
tion of the evolution of the specific star-formation rate with 
mass and redshift, and the evolution in the halo densities 
for a fixed Vcirc approximately cancel to predict that there 
should be very little evolution in the _B-band Tully-Fisher 
relation. 



3.1.2 The Stellar Mass Tully-Fisher Relation 

The rest-frame B-band Tully-Fisher relation is sensitive 
to recent star formation, and so perhaps the more fun- 
damental relation is the stellar mass Tully-Fisher relation 
(Mi versus Vasym) which refiects the time-integrated star- 
formation history. The 2 = stellar mass Tu l ly-Fisher rela- 
tion is well established (Bell fc do Jond l200ll : [p"izagno et al.l 
[2007; Mover et al. 2008), and in Fig. [6] we combine the stel- 
lar mass estimates with the inclination-corrected rotation 
speed of our galaxies to investigate how this varies as a 
function of redshift. Again, in this plot, we also include 
a number of hi gh-redshift measureme nts; at 2 ~0.6 an d 
2 ~ 1.3 from [Miller et all (|2011i ') an d (iMiller et al.ll2012l') 
the 2 = 1-3 lensed galaxy surveys ([Swinbank et al.l I2OO6I : 
iJones et al. I l2010bl : iRichard et aLlbOllI ) as well as the the 
2 ^2 SINS and z > 3 SIN FONI IFU surveys (|Cresci et all 
I2OO9I : iGnerucci et al.|[2oTll) . Where necessary, we have cor- 
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Fig ure 5. Left: The evolu tion of tlie rest-frame _B-band Tully-Fislier relation. We baseline our results against the z = measurements 
from lPierce k TuUvl l ll992l) . For the high-r edshift samples, we com bine the obser vations from this wor k with a number of othe r surveys. At 
z ~0 1-1 we include the field surveys from lBamford et al.l ll2005h (2 = 0.3 fi eld). IWeiner et al.l ||2006D (TKRS z = l fi eld) andlMiller et al.l 
ll201ll ') (z = 0.6 DEIMOS). At z ^ 1 we include the cluster arc survey s of [Sw inbank et al.l (|200^ ) (z = 1 arcs) and [joiies et al.l l l2010bl ') 
(z = 2— 3 arcs) and at z ^2 we include the SINS and AM AZE surveys iForste r Schrcibcr ot al. ( 20091 ) : [Cnerucci et al. I 1I2OIII) . The solid 
line denotes the z = B-band TuUy-Fisher relation from llPierce fc Tullvlll9 92). Right: The evolution of the zero-point of the B-band 
TuUy-Fisher relation. The open symbols denote individual galaxies and the solid points show the median (and error) in seven redshift 
bins betwe en 2 = and 2 = 3. The curves shows the pre dictions of the evolution of the _B-band TuUy-Fisher relation from semi-analytic 
models bv lBower et al ] l|2006l ) and iDutton et aP l|201ll ) which both predict evolution of AMb < 1 mags between 2 = and 2 = 3, less 
than seen in the observations. 



reeled Ihe samples to use V2.2 using V2.2/ Vmaj: = 0-93 if 
only V max is given in the literature (see also iDutton et al.l 
|2003)). 

In comparison to the rest-frame _B-band TuUy-Fisher 
relation, Fig. [5] shows that the the redshift evolution in the 
stellar mass TuUy-Fisher relation is weak. Due to the het- 
erogeneous selection functions of these various surveys, we 
again restrict the search to the evolution in the zero-point 
with redshift. We therefore construct the same seven redshift 
bins from z — 0-3 as in § 13.1.11 and measure the zero-point 
offset in each bin (using a non-linear least squares regression 
weighted by the velocity errors in each galaxy; (Fig. (H)). The 
data suggest that the zero-point offset from z = 2.5 to z = 
is AM*,z=o /Mj, = 2.0±0.4 (i.e. at a fixed circular velocity 
the stellar mass of galaxies has increased by a factor ~ 2.0 
between z — 2.5 and z = 0). The weak evolution in the stel- 
lar mass TuUy-Fisher relation (in particular below z ~ 1) is 
due to the fact that individual galaxies evolve roughly along 
the scaling relation rathe r than due to weak evolution in the 
galaxies themselve s (e.g. iPortinari fc Sommer- Larsenll2007l : 
iBrooks et al.|[201ll ). This indicates that the baryon conver- 
sion efficiency, rj— (M* /M200) / {^b / ^0, is fixed with red- 
shift. 

In Fig. |6] we also compare the zero-point evolution in 
the TuUy-Fisher r e lation to t hat from semi- a nalyt ic models 
from iBower et all (120061 ) and IDutton et al] l|201ll) and the 
cosmological hydrodynamic Galaxi es-Intergalactic M edium 
Interaction Calculation (GIMIC, ICrain et af] |2009| ). The 
iBower et al] (2006) and GIMIC simulations both use the 
same underlying dark matter distribution, although the 



GIMIC simulation follows the hydrodynamical evolution of 
gas, radiative cooling, star formation and chemo-dynamics, 
and includes a prescription for energetic feedback from su- 
pernovae. In contrast to simulations of small periodic vol- 
umes, GIMICs unique initial conditions enable it to fol- 
low, at high-resol ution, cosm ologicall y representative vo l- 
umes to 2 = (see lCrain et al . 2009; M cCarthy et al.ll2012l '). 
To search for evolution in the TuUy-Fisher zero-point in the 
simulations, we only consider galaxies with star-formation 
rates >1 Mq yr~^ (at all redshifts) so that a reasonable com- 
parison to observational data can be made. Both the semi- 
analytic and hydrodynamic simulations shown in Fig.[S]pre- 
dict strong evolution, with the GIMIC simulation suggesting 
that, at a fixed circular velocity the stellar mass of galax- 
ies should increase by 2.5 ± 0.5 x between z = 2 and z = 0. 
This evolution is m atched by the semi-analytic model from 
Dutton et al.l ( 201 j) . However, the semi-analytic model from 
Bower et al.l l|2006l ) predicts somewhat stronger evolution in 
the stellar mass for a fixed circular velocity. 



It is interesting to note that the GIMIC simulation pre- 
dicts that the scatter in the TuUy-Fisher relation should not 
evolve strongly with redshift (for a fixed circular velocity, 
the scatter is 25 ± 5% over the redshift range z = 0-3). Our 
dynamical me asurement s, particula rly when com bined with 
those from .Weiner et al.1 ( ,2006 ) and lMiller et al.l (2012), are 
sufficiently large in number to compare to the models, and 
Fig.[B]shows that the scatter in the observations is consistent 
with the models (for all but the last bin in redshift where 
the number of data-points is <8), , suggesting that the sam- 
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Figur e 6. Left: The evolution of the stellar mass TuUy-Fisher relation. We baseline the evolution against the z = work from lPizagno et alj 
1I2OO5I '). The high -redshift points are co mpiled from t he in termediate- and high- redshift (z ~0.6 DEIM OS and z ~1.3 DEIMOS) 
observations from | Miller et al. 1 1 I2OIII) and lMiller et al.lll2012D: the z = l and z = 2-3 cluster arc surveys from lSwinbank et al. 1 1I2OO6I') and 
Ijones et al] ll2010d) with s tellar masses from lRichard et al.l 1 I2OI1I) : and the z ~2-3.5 SINS and AMAZE surveys from lCresci et al.l 1I2OO9I') 
and [GneruccT'er'al] ll201lh . The symbols show individual galaxies. The solid line denotes the correlation at 2: = from iPizaenoet al.l 
ll200g ) (corrected to a Chabrier IMF). The dashed line is best-fitting zero point to the z = 2 sample galaxies (for a fixed V2.2) which 
shows an offset of of AM<,. z— p / M<: = 2 . ili 0. 4 between z = and 2 = 2.5. The do tted lines denote the TuUy-Fisher at 2 = 2 from the 
numerical simulations from lCrain et al.l l|2009j) and (see also iMcCarthv et al.ir2012h . which predict evolution in both the zero-point and 
slope. Here, we concentrate on the zero-point evolution, and note that the predicted evolution for a disk with circular velocity 100— 
200kms~^ is an increase in stellar mass of a factor 1.5-3 X between z = 2 and 2 = (equivalently, at high-redshift the maximum circular 
velocity is greater for the same stellar mass which may be consistent with the galaxies being compact at high-redshift and larger at 
low-redshift) . Right: The evolution of the zero-point of the TuUy-Fisher Relation. The symbols denote individual points (coded by the 
survey), whilst the solid symbols denote the average i n six redshif t bins. We also overlay the redshift evolution of t he zero-point of the 
TuU y-Fisher relat i on fro m the numerical model from ICrain et al. 1 1I2OO9I ) as well as the semi-analytic models from iBower et al.l 1I2OO6I ) 
and lDutton et akl 1I2OIII ) (DvdB). These galaxy formation models predict evolution in the zero-point of the TuUy-Fisher relation out to 
2 ~ 3 which is consistent with the observed trend given the large uncertainties in the latter. 



pies may now be sufficiently large (with rotation curves well 
enough sampled) that the scatter is intrinsic. 

3.2 The Redshift Evolution of the Mass-to-Light 
ratio 

Another route to examine the evolution in the B-band and 
stellar mass TuUy-Fisher relations is to combine the offsets 
and measure the evolution of the mean mass-to-light ratio 
with redshift. We caution that the conversion of zero-point 
offsets to offsets in mass-to-light assumes that star-forming 
galaxies form a homologous family and that the evolution 
is a manifestation of underlying relations between mass-to- 
light ratio and other parameters, such as star-formation his- 
tory, gas accretion and stellar feedback. 

Under these assumptions, in Fig. [7] we show the evolu- 
tion of the rest-frame B-band mass-to-light ratio. As ex- 
pected from the zero-point TuUy-Fisher offsets, this fig- 
ure shows that the average mass to light ratio of star- 
forming galaxies decreases strongly from z = up to 
z = 1 and then flattens. This behavior is consistent with 
the previous demonstration that star-forming galaxies at 
high-redshift have lower stellar masses and higher B- 
band luminosities. The strongest evolution occurs up to 
z ^1, AM* / Lb = 1.1 ± 0.2, consistent with previous studies 



(|Miller et al.ll201ll ) (the fractional change in mass-to-light 
ratio over this period is A(M /Lb) / (M /Lb)z=o ~3.5 be- 
tween 2 = 1.5 and z = 0, with most of the evolution occuring 
below z — 1). 

We note that we examined whether the evolution in 
the mass-to-light ratio could be reproduced using simple 
star-formation histories, ranging from i) a constant star- 
formation rate (with a formation redshift, Zf = 4-8); or ii) a 
set of exponentially decreasing star-formation rates with e- 
folding times ranging from 0.25-10 Gyr (and formation red- 
shifts ranging from Zf = 2-10). However, using these simple 
star-formation histories (adopting a Chabrier IMF with 0.5- 
1 solar metallicities and the Padova (1994) stellar evolution 
tracks), we are unable to find a acceptable fit with a sin- 
gle star-formation history. This could be because the "av- 
erage" star-formation history is more complex than a sim- 
ple star-formation model, or because the current low and 
high-redshift data can not be linked by a simple evolution- 
ary model. However, in Fig. [7] we also overlay the predicted 
evolution of the B- band mass-to -l ight r atio f rom the semi- 
analyt ic models of iBower et al.l (|2006l ) and iDutton et al.l 
(j201ll ) , both of which predict a sharp decline to 2 ~ 1 and 
then flattening to higher redshift, which provides a reason- 
able match to the observations. 
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Figure 7. The B-band mass-to-light ratio as a function of red- 
shift for star-forming galaxies from this work plus the compila- 
tion of studies shown in Fig. \5\ and \E\ The dashed- and dot- 
dashed- lines show the predicted evolu tion of the mass-to -light 
ratio as a function of redshift from the iBower et al. I ll2006l ') and 
iDutton et al ] l|201ll l semi-analytic models. This figure shows that 
in both the observations and galaxy formation models there is 
strong evolution in the mass-to-light ratio of galaxies with red- 
shift. The strongest evolution in the data is seen for z = 0-1, with 
significant fiattcning above z = l. 



3.3 The Evolution of the Size— Rotational Velocity 
Relation 



In a ACDM cosmology, the sizes of galaxy disks and their 
rotational velocities should be proportional to their parent 
dark halos, and since the halos are denser at higher redshift, 
for fixed circular v elocity the disk sizes should scale inversely 
with Hubble time ||Mo et al.lll998l ). In this scenario, disks at 
z = l and z = 2 should be 1.8 and 3.0 x smaller respectively 
at fixed circular velocity. With the high-resolution measure- 
ments for the galaxies in our sample we can measure the 
extent of the Ha (and together with the dynamics) compare 
these relations at 2 1.5 to star- forming disks at z — 0. 

In Fig.|8]we show the relation between the Ha half-light 
radius and rotational velocity for the galaxies in our sample 
(using the asymptotic velocity as measured from the disk 
modeling). In this plot, we also include the me asurements 
from the SINS galaxies from lCresci et al] (|2009l ). First, we 
note that the Ha half-light radii for the SHiZELS galax- 
ies {rh = 2.4 ± 0.8 kpc) arc comparable to the average rest- 
frame UV/optical h alf light radii of star - forming galaxies 
at this redshift (e.g. [Conselice et al.ll2003l : iLaw et al.ll2007l : 



the errors on each data-point) suggests that this offset is 
a factor Arh / rh,z=o = 1-4 ± 0.3, which is consistent (given 
the smaller number statistics in the high-redshift measure- 
ments) with that pre dicted fr om halo-disk scali ng argu- 
ments ||Mo et al.lll998l : see also button et al.ll201ll ). 

The product of the rotational velocity and size of the 
disk define the specific angular momentum, and the evolu- 
tion (at a fixed circular velocity) with redshift also encodes 
^.0^ information on the relation between the star-formation 
processes in the disk and the halo spin par ameters (e.g. 
iNavarro fc Steinmetj [20001 : iTonini et al.ll2006l ). The offsets 
seen in the ri/2-Vasym relation compared to local scaling 
relations in Fig. [8] translate directly into offsets in the 
specific angular momentum for disks between z = and 
2; ~2. Indeed, the offsets in ri/2-Vasym suggest that, for 
fixed circular velocity the specific angular momentum of 
disks at z = 2 is a factor 1.4 ±0.3 lower than those at z = 
(INavarro fc Steinmetj I2OO0D . Thus, if high-redshift disks 
evolve into local disk galaxies, then the specific angular 
momentum must be increased, either by removing material 
with low specific angular momentum through outflows, or 
the galaxy disk must decouple from the halo and acquire 
angular momentum from accreting gas. 

Finally, in Fig. [S] we also show the relation between 
stellar mass and Ha half light radius for the g alaxies in our 
sample compared to the z = relation from iDutton et al.l 
(2007). The high-redshift galaxies are offset from the local 
relation by Arrh / rh = 1-2 ± 0.3 for a fixed stellar mass. Al- 
though this is consistent with no evolution, we note that the 
predicted evolution for a fixed stella r mass between z = -2 
from the semi-analytic model from IDutton et al.l (|201ll ) is 
Ar,h/rh = 1.2. 



ISwinbank et al. I I2OI0I ) 



In order to compare to z = data, w e show the mea- 
surements of local disks from INavarro fc Steinmctz (2000) 
(whic h is based on a com pilation from iMathewson et al.l 
Il992l and ICourteaul Il997l ). but converting the disk scale 
length (rd) to half-light radius using ri/2 = l-68rd, as ap- 
propriate for a disk with an exponential profile. As can be 
seen from this figure, the high-redshift galaxies have half- 
light sizes (deconvolved for seeing) which are smaller for 
fixed circular velocity than those at 2 = 0. Fitting the lo- 
cal the 2 = relation to the high-redshift data (only al- 
lowing a zero-point shift and inverse weighting the fit by 



3.4 Metallicities and Spatially Resolved Chemical 
Abundances 

The internal enrichment (and radial abundance gradients) 
of high-redshift star-forming galaxies provides a tool for 
studying the gas accretion and mass assembly process such 
as gas exchange (infiows/outflows) with the intergalactic 
medium. However, the sm all sizes of high redshift galaxies 
(typi c ally rh ~ 2-3 kpc; e .g. Law et al. 20071 : Conselice et al.l 
|2008l : ISmail et al. I I2OO7I : ISwinbank eTsI} I2OIOI ) mean that 
sub-kpc resolution observations are required to measure the 
abundances reliably. At high redshifts the gas-phase metal- 
licity of individual galaxies can only be measured with lim- 
ited accuracy, usually through the [Nil] / Ha emission-line 
ratio as an indicator of oxygen abundance and calibra- 
tions between gas phase emission line an d stellar absorp- 
tion m etallicities, such as those obtained bv lPettini fc Pagell 
l|2004l) . Recently, the first reliable spatially resolved mea- 
surements of [Nil] / Ha have been made in two highly am- 
plified galaxies at 2 ~ 1.5-2.5, measuring strong, negative 
gradients and suggesting that if 2 ~ 2 starbursts evolve 
into spiral disks today, then the gradients must fiatten by 
A log(0 / H) / AR~ 0.05-0.1 dex kpc over the last 8-10 Gyr 
l|jones et al.l2010al : lYuan et al.ll201ll ). Other studies at lower 
resolution have derived weaker gradients that are negative, 
or consistent with z ero for rotating and interacting galaxies 
(|Cresci et al.ll2010l : lOuevrel et al.] |2012l ). Positive gradients 
for interacting systems may be a signature of the redistribu- 
tion of the metal-rich gas produced from a central starburst; 
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Figure 8. Left: The relation between the circular-velocity and Ha 
half-light radius for high-redshift galaxies. We include the results 
from t he SHiZELS galaxies in our sample as well as those from 
SINS llCresci et al. 2009). T he so lid line shows the 2 = relation 
from [Navarro fc Steinmet3 ll2000l ') ( which is based o n a compila- 
tion from lMathewson et al.lll992l and lCourteaulll997h . Fitting the 
local relation (assuming a fixed slope) to the high-redshift points 
suggests and offset Arjj / z— g = 1.4 ± 0.3, which is shown by the 
dot-dashed line (the scatter is shown by the dotted lines). This 
evolution is c onsistent with that p redicted by halo— disk scaling ar- 
guments (e.g. lDutton et al]|201ll : blue dashed line). Right: the re- 
lation between stellar-mass and half-light radius in the SHiZELS 
and SINS surveys. The solid line shows the z = relation from 
iDutton et al. 

1 ll201ll ) whilst the dot-dashed curve shows the offset 
of the high-redshift measurements from the local relation (as- 
suming a fixed slope), which shows a Ar^ii / rji = 1.2 ± 0.3 offset 
for a fixed stellar mass. We also show the predic ted evolution 
oiz = 2 galaxies from the 2 = relation from the iDutton et al.l 
ll201lh semi-analytic model. 
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l|Werk et al.ll2010l : lRupke et al.ll2010l ). or possibly as a result 
of inflow of relatively unenriched gas from the halo (or IGM) 
if it is able to intersect the central regions of the disk galaxy 
without being disrupted. 

With our high-resolution (and high signal-to-noise) 
resolved observations we can measure the spatially re- 
solved emission line ratios and hence the distribution of 
chemical abundances. First, we note that the average 
galaxy integrated [Nil] / Ha ratio for the galaxies in our 
sample is [Nil] / Ha = 0.27 ± 0.09, which corresponds to 
12 + log(0/H) = 8.58 ±0.07 (equivalently, ^ 0.73 ± 0.16 Zp 
wher e Zq is 12 + log(0 / H) = 8.66 ± 0.05; lAsplund et all 
12004 ). We caution that the [Nil] /Ha ratio is sensitive 
to shock excitation, the ionization state of the gas, and 
the hard ionization radiation field from an active galac- 



tic nucleus (iKewlev fc Dopital |2002| : iKewlev et all I2OO6I : 
iRich et al.ll201Gl ). However, our galaxy- integrated [Nil] / Ha 
ratios are consistent with star formation throughout the disk 
(Fig.©. 

Of the nine galaxies in our high-redshift sample, 
seven have [Nll]A6583 emission which is sufficiently strong 
(S/N>10 in the combined spectra) to allow spatially re- 
solved studies of the chemical abundance gradients. First 
we subtract the large-scale velocity motions from the data- 
cube (as defined by the best-fit dynamical model). We then 
use the disk inclination and position angle to define ~lkpc 
annuli (centered at the dynamical center) and integrate 
the spectrum in each annulus to measure the [Nil] / Ha ra- 
tio. In Fig. [9] we show the Ha emission line maps as well 
as the spatially resolved [Nil] / Ha ratio and overlay the 
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Figure 9. Metallicity gradients for the seven galaxies in our sam- 
ple where spatially resolved measurements can be made. Left: Ve- 
locity integrated Hct emission line maps. The dashed annuli de- 
note the ~ 1 kpc regions from which the [Nii]/Hq were extracted. 
Right: The [Nil] A 6583 / Ha emission line ratio as a function of 
physical radius. We overlay the best-fit linear regression and we 
report the gradients for individual galaxies in Table 1. In all cases 
the gradients arc negative or consistent with zero, as expected for 
inside-out growth. 
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Figure 10. Metallicity gradients versus redshift for the galax- 
ies in our sample compared to the Millty-Way and other low- 
and high-redshift observatio ns. The Milky-Way t h in and thic k 
disk measurements are fromlGilmore et al.l lll995^ : iBelll l ll996t) : 
iRobin et all lll996l '): lEdv ardsson et al.! ('l99^); the low-redshift 
star- forming gal axies are fr om Ruokc ct al. ( 20101^ and th e z=l - 
1.5 data from lYuan et al.l ||2011| ') and iQuevrel et al.l | |2012|) . 
The median abundance gradient for our SHiZELS sample is 
Alog(0 /H) / AR=-0.027±0.005dexkpc-i which is consistent 
with the thick disk of the Milky- Way. We also show the theoretical 
evolution of the metallicity gradient with redshift from the GIMIC 
simulation llCrain et al. 1 120091 : iMcCarthv et al1l2012f) . The dark 
grey shows the range of metallicity gradients for all disk galax- 
ies in the simulation in the mass range 9.5 < log(M / M©) < 11.5, 
whilst the light grey denotes the Icr scatter at each epoch. The 
GIMIC simulation predicts the gradients should be steep around 
2 ~3 where gas accretion from the IGM is at its most efficient, 
but should flatten by a factor ~2 X between 2 = 3 and 2 = 0. 



best-fit linear regression in eacli case. We report the gra- 
dients individually in Table 2. In all cases, the metallic- 
ity gradients are negative or consistent with zero, with an 
average A log(0 / H) / AR= -0.027 ± 0.005 dexkpc"\ and 
similar to those recently foun d for a sample of z ~ 1 
starbursts l|Quevrel et al.ll2012l ). This average gradient is 
steeper, but consistent with the thick disk of t he Milky- Way, 
Alog(O/ H)/AR=-0.01±0.02dexkpc~^ (iGilmore et al.l 
11995 



Beii|[l99^ : lRobin et al. I ll996l : lEdvardsson et al.lll993h . 

which has a mass weighted age of 8-12 Gyr (equivalently, a 
form ation redshift z = 2-4; iBensbv et al.ll2007l : iRuchti et al] 

In Fig. [To] we compare the metallicity gradients to those 
from the GIMIC simulation (which include cold gas accre- 
tion from the IGM). We calculate the disc phase metallicity 
in cylindrical coordinates (along the disc) at z — 0, 1, 2 & 
3 in radial bins of width 1 kpc for all galaxies in the sim- 
ulation with masses between M* = lO^ '^"" Mq. At z ~ 2, 
the model disk galaxies from the GIMIC simulation display 
negative metallicity gradients which are consistent with our 
observations. Moreover, between 2 = 2 and z = the GIMIC 
simulation predicts that the gradient should flatten by a 
factor ~ 2. In the simulation, this is a consequence of the 
outer disk enrichment increasing with decreasing redshift 
(as the gas inflow rates decline and gas is redistributed) 
progressively flattening the radial gradient within the thick 



di sk. This is qualitat ively similar to the numerical models 
of iDekel et al.l l|2009al ) and where the "cold streams" (which 
are most prevalent around z ~2) penetrate the halo and 
intersect and deposit rel atively unenriched gas at onto the 
galaxy disks at ~15kpc (jPekel et al.ll2009al lbl). 



4 CONCLUSIONS 

We have presented AO-assisted, spatially resolved spec- 
troscopy of nine star-forming galaxies at z = 0.84-2.23 se- 
lected f rom the UKIRT, w ide-field narrow-band, HiZELS 
survey (|Sobral et al.ll2012bh . Our main results can be sum- 
marised as: 

• Using Ha emission line dynamics, we find that the 
ratio of dynamical-to-dispersion support for the sample is 
'^max sin(i) / (T = 0.3-3, with a median of 1.1 ±0.3, which 
is consistent with similar measurements for both AO and 
non-AO studies of star-form ing galaxies at this epoch(e.g. 
iForster Schreiber et al.ll2009l ). At least six galaxies have dy- 
namics that suggest that their ionised gas is in a large, rotat- 
ing disk (in at least four of these we detect the turn over in 
their rotation curves). For these galaxies that resemble ro- 
tating systems, we model the galaxy velocity field and derive 
the inclination-corrected, asymptotic velocity. 

• We combine our dynamical observations with 
a number of previous studies of inte rmediate- and 
high-redshift galaxies (in particular from Bamford et alJ 
2005; Woincr ct al. 2006; Forstcr Schreiber ct al. 2008 



Jones et ani2010bl : iMiUer et al.ll201li and iMiUer et al.ll20r3 ) 



to investigate the evolution of the TuUy-Fisher relation. We 
show that there is strong evolution in the zero-point of the 
rest-frame B-band TuUy-Fisher relations with redshift such 
that at a fixed circular velocity, the B-band luminosity in- 
creases by A Mb = 2 mags up to z ^ 2 (a change in lumi- 
nosity of a factor ~ 6; e.g. see also iBamford et al. I l2005l : 
IWeiner et allbOOel : iMiller et al.ll201ll ). 

• In contrast, the stellar mass TuUy-Fisher relation 
shows much more modest evolution; for a fixed circular ve- 
locity the average stellar masses has increase d by a factor 
2.0± .4 between z = 2 a nd z = (e.g. see also ICresci et al.] 
I2OO9I : iMiller et al1l2012l ). The weak evolution in the stel- 
lar mass TuUy-Fisher relation (in particular below z r^l) 
is due to the fact that individual galaxies evolve along the 
scaling relations rather t han due to weak evolution in the 
galaxies themselve s (e.g. iPortinari fc Sommer-Larsenll2007l : 
iBrooks et aLll201ll : iDutton et al.ll201lh . and indicates that 
the baryon conversion efficiency, 77 = (M* / M200) / (f^b / ^^o), 
is fixed with redshift. The observations also suggest that the 
scatter in high-redshift TuUy-Fisher relation is comparable 
to that predicted by the models, suggesting that the sam- 
ples may now be sufficiently large (with rotation curves well 
enough sampled) that the scatter is intrinsic. 

• Combining the rest-frame B-band and stellar 
mass TuUy-Fisher relations we show that the mass-to- 
light ratio of star-forming galaxies evolves strong with red- 
shift, decreasing by AM/Lb = 1.1±0.2 between z = and 
z=l, but then ceases to evolve further above z = l. This 
change in mass-to-light ratio in the B-band is a factor 
A (M / Lb) / (M / Lb)z=o ~ 3.5. We show that the evolution 
in the TuUy-Fisher relations and mass-to-light ratio is in line 
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with both cosmologically-bascd hydrodynamic simulations 
and semi-analytic models of galaxy formation. 

• Using the Ha half light radius and circular veloc- 
ity, we also find that the SHiZELS galaxies are a factor 
1.5 ± 0.3 X smaller for fixed circular velocity than for disks at 
z = 0. If the z ^2 and z = disks are linked in a simple evo- 
lutionary scenario then high-redshift disks must acquire spe- 
cific angular momentum, either by removing material with 
low specific angular momentum (e.g. through outflows), or 
through acquisition of angular momentum from accreted gas 
at late times. 

• Finally, we measure the spatially resolved [Nil] / Ha 

to measure the mctallicity gradients. In all cases, the motal- 
licity gradients are negative or consistent with zero, with 
an average Alog(0 / H) / AR= -0.027 ± 0.005 dexkpc"\ 
which is consistent with that seen in the thick disk of the 
Milky- Way. We show that these metallicity gradients are 
comparable to those predicted for the gas disks of star- 
forming galaxies in the cosmologically based hydrodynamic 
simulations. In these models, the inner disk undergoes initial 
rapid collapse and stax-formation with gas accretion (either 
from the halo and/or IGM) depositing relatively unenriched 
material at the outer disk, causing negative abundance gra- 
dients. At lower redshift (when gas accretion from the IGM 
is less efficient) the abundance gradients flatten as the outer 
disk becomes enriched by star- formation (and/or the redis- 
tribution of gas from the inner disk). 

Overall, wc demonstrate that well-sampled, dynami- 
cal measurements of high-redshift star-forming galaxies can 
constrain the rotation curves, measure the evolution of the 
basic disk scaling relations and their radial abundance grar 
dients. We show that the current observational samples with 
well resolved data arc now approaching sufficiently large 
samples that it is possible to refine or refute models of disk 
galaxy formation. 
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